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Abstract 

We examine the maximum entanglement reached by an initially fully aligned state evolving in 
an XY Heisenberg spin chain placed in a uniform transverse magnetic field. Both the global 
entanglement between one qubit and the rest of the chain and the pairwise entanglement between 
adjacent qubits is analyzed. It is shown that in both cases the maximum is not a monotonous 
decreasing function of the aligning field, exhibiting instead a resonant behavior for low anisotropics, 
with pronounced peaks (a total of [n/2] peaks in the global entanglement for an n-spin chain), 
whose width is proportional to the anisotropy and whose height remains finite in the limit of small 
anisotropy. It is also seen that the maximum pairwise entanglement is not a smooth function of the 
field even in small finite chains, where it may exhibit narrow peaks above strict plateaus. Explicit 
analytical results for small chains, as well as general exact results for finite n-spin chains obtained 
through the Jordan- Wigner mapping, are discussed. 

PACS numbers: 03.67.Mn, 03.65.Ud, 75.10.Jm 
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I. INTRODUCTION 



Quantum entanglement has been long recognized as one of the most fundamental and 
intriguing features of quantum mechanics l|. It denotes the ability of composite quantum 
systems to develop correlations which have no classical counterpart. Interest on all aspects 
of entanglement has grown enormously since its potential for permitting radically new forms 
of information transmission and processing was unveiled 
an essential resource in the field of quantum information science 



glement measures have been introduced 



5], being now considered 
6], where rigorous entan- 



7J, |8| . The interest has also extended to other areas 



like condensed matter physics, where it has provided a novel perspective for the analysis of 

12,ll3|. 

Spin systems with Heisenberg interactions [14j, ll5j constitute a particularly attractive 
scenario for studying quantum entanglement. They provide a scalable qubit representation 



suitable for quantum processing tasks 



;iement. iney 

3,h, fly 



and can be realized by diverse physical 



systems such as cold atoms in optical lattices [20|], quantum dots [la, ll7[ and Josephson 



junctions arrays 



211 ] . Accordingly, several investigations of entanglement in ground and 



thermal equilibrium states of Heisenber g sp in chains subject to an external magnetic field 



have been made (see for instance 



0,Q, 



ii 



12 



2i,|23|,|2j). There 



rave also been relevant 



25 



26 



28 



29|), which 



studies of entanglement dynamics in spin chains (for instance 
discuss in particular the evolution of initial Bell states and the ensuing "entanglement waves" 
251 ] . non-ergodicity and dynamical phase transitions starting with equilibrium states 26], 



decoherence waves 



271 ]. evolution in varying magnetic fields 28j], generation of cluster states 



29] as well as other issues. 

In the present work we want to focus on a particular aspect, namely the generation of 
entanglement in an interacting spin chain with fixed parameters starting from an initially 
fully separable aligned state, and examine the maximum entanglement that can be reached 
as a function of the anisotropy and the uniform transverse magnetic field (control parame- 
ter). We will concentrate here on the glo bal entanglement between one qubit and the rest of 
the chain and on the pairwise entanglement between neighboring qubits, within the context 
of a cyclic XY chain with nearest neighbor interactions [14]. Questions which immediately 
arise include the possible existence of a threshold anisotropy for reaching maximum global 
entanglement (saturation), the maximum pairwise entanglement that can be reached and, 
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most important, the behavior with the applied magnetic field. It will be shown that con- 
trary to what can be naively expected, the maximum global entanglement reached is not a 
monotonous function of the aligning field, but exhibits instead a typical resonant behavior for 
low anisotropies, with narrow peaks located at characteristic field values, entailing a high 
sensitivity suitable for entanglement control. The pairwise entanglement exhibits a more 
complex resonant response, since it is affected by a competition between two incompatible 
types (essentially of positive or negative spin parity). These resonances remain finite in the 
limit of vanishing (but non-zero) anisotropy in finite chains, considering sufficiently long 
time evolutions. On the other hand, for large anisotropies they merge into a single broad 
maximum centered at zero field, with global saturation reached within a field window. 

Our results are based on a fully exact treatment of the finite n-spin chain based on the 
Jordan- Wigner transformation [141 ]. explicitly verified for the case of two and three-qubit 
chains. The Hamiltonian and the entanglement measures employed are discussed in section 
II. Section III contains the results, discussing first the two and three-qubit cases and then 
the exact results for general n-qubit chains. Finally, conclusions are drawn in IV. 



II. FORMALISM 



We consider n qubits or spins in a cyclic chain interacting thr oug h an XY nearest neighbor 
coupling, embedded in a uniform transverse magnetic field {14I Iis| . The Hamiltonian reads 

n 

H = bS z — + v yS y j+1 ) (la) 

3=1 

1 n 

= hSZ ~2 X>4 s 7+i + 944+1 + h - c -) » ( lb ) 

3=1 

where S z = YTj=i sZ j ls the total spin along the direction of the magnetic field b, v,g = 
(v x ±v y )/2 and n + 1 = 1. We will consider the evolution of the state which is initially fully 
aligned antiparallel to the magnetic field, 

Mt)) = eM-m\i ■ ■ ■ i) , (2) 

where t denotes time over h, and examine the emerging global entanglement between one 
qubit and the rest of the chain, as well as the pairwise entanglement between contiguous 
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qubits, arising for non-zero anisotropy 7 = g/v (for g = the initial state is an eigenstate 
of H and hence no entanglement is generated). 

Q 

Since we are dealing with a pure state, the first one is determined by the entropy [7[ 

Ex = -Tr px log 2 pi , (3) 

of the reduced one-qubit density p\ = Tr n _xp, where p = \ty (t)\ is the full density 
matrix, with E\ = for p\ pure (p\ = p\) and E\ — 1 (maximum) for pi fully mixed. The 
second one is the entanglement of formation [7| of the adjacent pair density p 2 = Tr n _ 2 p, 
which can be calculated as 

E2 = -^2<b> !og 2 Qv , (4) 

u=± 



where g ± = (1 ± y/l - CD/ 2 and 

C 2 = Max[2A m -Tr#,0], R = ^f 2 , (5) 



is the concurrence [8j, with A m the greatest eigenvalue of i? and p 2 = 4sJs^ +1 p 2 s^ +1 sJ the 
spin-flipped density. It satisfies < Ci < 1. Since tracing out qubits of a subsystem can 
be considered a LOCC (local operations and classical communication) transformation, it 
cannot increase entanglement [7| and hence E 2 < £/i, with E7 2 = for a pure two qubit 
state (in which case q± become the eigenvalues of pi). 

As E 2 is just an increasing function of C 2 , pairwise entanglement is usually directly mea- 
sured through the latter, which is more suitable for analytic description. The corresponding 
measure of the global E\ entanglement is the square root of the one-tangle {sj], 

Ci =2VDetpi = ^(l-Trp 2 ), (6) 

which coincides with C 2 for a pure two qubit state and satisfies G\ > C 2 in the general case 

with Cij the concurrence of the 



(actually the more general inequality C{ > ^Ylj^i ^ 
pair and Cf the one-tangle of qubit i, conjectured in [30(], was recently proven [3l|]). Both Ei 
and C\ are measures of the disorder associated with p\ and are hence increasing functions 
of one another. 

Due to the symmetries of H and the present initial state, |^/(t)) will be invariant under 
translation (j — > j + 1) and inversion (J — > n + 1 — j), and will have positive spin parity 
P = exp[in(S z + n/2)], as this quantity is preserved by H ([H,P] = 0). The reduced 
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density ps = Tr n _sp of any subsystem S will then depend just on the distance between its 
components and will commute with the subsystem parity Ps = IX/eS ex P[* 7r ( S 2 + 1/2)], as 
the reduction involves just diagonal elements in the rest of the chain. In the case of pi, this 
implies that it will be the same for all qubits and diagonal in the standard basis |l) of 
s z eigenstates: 

'pit) 
l-p(t) 

where p(t) represents the one-qubit spin flip probability 



Pi 



(7) 



p(t) = {s*) t + 1/2 =(S') t /n + 1/2, (8) 
(here (0) t = (t)\0\^ (t)) and spin operators are considered dimensionless) . Hence, 

C 1 (t)=2y/p(t)[l-p(t)], (9) 

with d(t) = 1 when p(t) = 1/2. 

The same symmetries lead to a pair density of the form 

f Px(t) a*{t)\ 



P2 



(10) 



p 2 {t) (3{t) 
p(t) p 2 {t) 
\a{t) p 3 (t) J 

in the standard basis |TT), III), lit), |U), where p 1 (t) + 2p 2 (t) +p 3 (t) = l, Pi(t) +p 2 (t) = p(t) 
and 



Pl (t) = ((4 + l/2)(^ +1 + l/2)) t , 
for adjacent qubits. Eq. (jSJ) becomes then 



C 2 (t) = 2Max[|a(*)| -p 2 (t), \/3(t)\ - y/ P i{t)p 3 {t),0] 



(11a) 
(lib) 

(12) 



where only one of the entries can be positive (this follows from the positivity of p 2 , which 



requires \a(t)\ < Pi{t)p 3 {t) , \/3(t)\ < p 2 {t)). Two kinds of pairwise entanglement can 
therefore arise: type I (|a(£)| > p%{t)) and type II (|/3(t)| > Pi{t)p${t)) , which cannot 
coexist and can then be present just at different times, and which stem from the positive (I) 
and negative (II) parity sectors of p 2 . 



The eigenvalues of H and the entanglement of its eigenstates are obviously independent 
of the sign of g, and for even chains also of the sign of v, as for n even it can be changed by a 
transformation s*' y — > (— l) J 'sJ' y . Due to time reversal symmetry, the emerging entanglement 
in even chains will then be also independent of the sign of b, while in odd chains that for 
(—6,1;) will coincide with that for (6, —v). We will then set in what follows v > 0, g > 0, 
and consider both signs of b. 

III. RESULTS 

A. Two qubit case 

Let us first analyze this simple situation, which nonetheless provides already some insight 
on the behavior for general n. Here C\ — Vi. The evolution subspace is spanned by 
the states |||), |||), and the pertinent eigenstates of H are |±) = m t |||) =F ^±|TT), with 
energies E± = ±A, where u± = y/(\±b)/(2\) and A = ^b 2 + g 2 . The state (J2J will then 
be independent of v and given by 

!*(*)> = E^^IIDI^) 

u=± 

= (cos At + if sin At) III) + if sinAtltt), (13) 

so that the spin-flip probability p(t) is 

p(t) = ¥ ^— 2 sin 2 \t. (14) 

Its maximum p m = g 2 / (b 2 + g 2 ) is thus a Lorentzian of width g centered at b = 0, satisfying 
Pm > 1/2 if |6| < g. Hence, for any g > the system will always reach maximum entangle- 
ment C\ = 1 within the field window |6| < g, at times t m such that p(t m ) = 1/2, where Eq. 
( TT3l) becomes a type I Bell state: 

|*(t ro )> = ±i(|TT) + e^\U))/V2 , cos0 = b/g . 

The maximum concurrence reached (Fig. [1]) is then 

c "' = c ™ = h;, ! s 'f! ■ d5) 

I s 2 + 1 , |S| ^ i 
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which is higher than the concurrence C = g/X of the Hamiltonian eigenstates V b ^ 0, 
becoming « 2g/\b\ for |6| ^> g. Ci{t) will follow the evolution of p(t) if p m < 1/2 (|s| < 1), 
but will develop saturated maxima plus an intermediate minima when p m > 1/2. 

We also note that for 6 = 0, i.e., where the gap E + — E__ = 2A is minimum and vanishes 
for g — > 0, maximum entanglement can in principle be attained for any g > 0. In this case 
the eigenstates |±) become independent of g and maximally entangled, and none of them 
approaches the aligned initial state for g — > (in contrast with the behavior for 6 7^ 0). 
The initial state becomes then equally distributed over both eigenstates (u± = 1/V2) V 
g > 0, implying |\?(t)) = cosgt|||) + isingt\ ff). Hence, in this case the only limit for 
reaching maximum entanglement (sin 2 gt = 1/2) for arbitrarily small but non-zero g is the 
long waiting time (t m = n/(Ag)). We will see that an analogous situation will occur for any 
n at particular field values. 



B. Three qubit case 

For n = 3, the evolution subspace is still two-dimensional and spanned by \ III) and 
the W-state M \W) = ( | X T T> + I TIT) + |tti»A/3, which for g = have energies -36/2 
and 6/2 — v. The coupling induced by g leads to eigenstates |±) = u T \ HI) =F u±|W) 
with energies E± — e ± A, where u± = a/[A ± (6 — v/2)]/ (2A), e = —(6 + v)/2 and A = 



y/(b-v/2) 2 + 3g 2 /4. We then obtain 

= e - fet [(cos At + sin At) HID + sin Xt\W)} 

which leads to 

pit) = —rr, rrz t"tt sin 2 At . (16) 

Fy ! 2[(6- v/2) 2 + 3^ 2 /4] v ; 

Its maximum p m = g 2 / (2A 2 ) is again a Lorentzian of width proportional to g but centered at 
6 = v/2 due to the hopping term, where p m = 2/3 (the value at the W-state), with p m > 1/2 
for |6 — v/2\ > g/2. Hence, for any g ^ there is again a field interval where maximum E\ 
entanglement is attained. The maximum of Ci(t) (Fig. [TJ top right) is then 

fl, \s\ < 1/2 

C ™=\ v^Ti , l>1/9 ," = {b-v/2)/g. (17) 

For |6| ^> v,g, C™ ~ V^gr/|6|, an asymptotic result which turns out to be valid V n > 3. 
The evolution of Ci(t) remains qualitatively similar to that for n — 2. Note also that for 



b/g (b-v/2)/g 




FIG. 1: (Color online). Top: Left: Maximum entanglement (measured by the concurrence) reached 
by the two qubit chain as a function of the (scaled) magnetic field for an initially aligned state. 
The dotted line depicts the concurrence of the Hamiltonian eigenstates. Right: Maximum global 
concurrence Cf 1 between one-qubit and the rest (upper curve, in blue) and maximum pairwise 
concurrence C™ (lower curve, in red+dashed pink) in the three qubit system, in terms of the 
(shifted+scaled) magnetic field. exhibits a sharp type II resonance at b = v/2. Bottom: Left: 
Plot of C\ and C 2 in the three qubit chain in terms of the spin flip probability p (0 < p < 2/3). 
Right: The temporal evolution of C\ and C 2 in the three qubit chain at the C™ plateau (6 = 
v/2 ± 0.65, left) an d at resonance (b = v/2, right). T = 2ir/\ is the period. Type I (II) sectors in 
C 2 are depicted in solid red (dashed pink) lines. 

b = v/2, i.e., where the gap 2A is minimum and vanishes for g = 0, maximum C\ is again 
reached for any g > 0, the situation being similar to that for n = 2 at b = 0. 

The behavior of the pairwise entanglement is, however, more complex. The W-state 
contains type II pairwise entanglement, but \^{t)) will first develop that of type I, so that 
transitions between both types can be expected to occur in the evolution for large g. From 
the expression of we obtain |a(t)| = y/p(t)(2 - 3p(t))/2, p 2 (t) = pi(t) = (3(t) = 

p{t)/2, so that Eq. ( fT2l) becomes 

C 2 (t) = \y/p(t)[2-3p(t)]-p(t)\ , (18) 
8 



which corresponds to type I (II) for p(t) < 1/2 (> 1/2). It is thus a non-monotonous 
function of p = pit) (left bottom panel in Fig. [1]), having a maximum at p — 1/6 (where 
C*2 = 1/3), vanishing at the "critical" value p = 1/2 (where C\ is maximum) and increasing 
again for p > 1/2 up to its absolute maximum at the endpoint p = 2/3, where Ci = 2/3 (i.e., 
the value at the W-state). Hence, saturation (C2 = 1) cannot be reached. Moreover, it is 
verified that C 2 {t)/C 1 {t) < l/y/2 (the maximum ratio allowed by the generalized inequality 
[30I for C12 = C13), the maximum reached for p —>■ or p ^ 2/3. The evolution of (^(i) 
will then not follow that of p(t) or C\(t) if p m > 1/6, developing for p m < 1/2 a minimum 
when p{t) is maximum, which will evolve into two vanishing points plus a type II maximum 
if p m > 1/2 (see right bottom panel in Fig. [T]). The maximum of C%{t) is then 



' 1/2-H 

s 2 +3/4 



\s\ < S c 

b-v/2 

1/3 , s c < \s\ < 3/2 ,8= — (19) 



s 2 +3/4 



\s\ > 3/2 



where s c = \/3 — 3/2 w 0.23 determines the second point where C 2 = 1/3 and encloses the 
region of dominant type II entanglement. It then exhibits a sharp type II peak at b = v/2, 
above a strict type I plateau (see Fig. [2]). Note that at b = v/2, C™ = 2/3 for any g > 0, as 
in this case the system will always reach the W-state if the waiting time is sufficiently long 
(t m = 7r/(v^3fi0)- For \b\ ^> v,g, C™ ~ g/\b\ ~ C™/V2, an asymptotic result which is again 
valid Vn>3. 



C. General n 

By means of the Jordan- Wigner transformation we may exactly convert the Hamil- 
tonian (TjQ) within a fixed spin parity subspace (P = ±1) to a quadratic form in fermion 
operators cL Cj, defined by cj = exp[— in J2i=i s i~ s 7)- ^ or a fi n he cyclic chain with posi- 
tive parity P — 1, the result for H' = H + bn/2 is 

i=i 

= J^(6 - u cos w fc )cj t4 - |^ sin u k i4c'l k + c!_ fc t4) , (20b) 
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where the fermion operators c' k , c k are related to Cj, cj by a finite Fourier transform 

with k half-integer for the present cyclic conditions: = — -^^p-,..., 1 ^ for n even and 
k = — ~ + 1, . . . , ~ for n odd. We then obtain the diagonal form 

H' = ^ A fc 4 a fc - | [*k - {b - v cos u k )) , 



Afe = y (6 - wcosu; fc ) 2 + 5' 2 sin 2 cu fc , (21) 

by a means of a BCS-like transformation c k = Ukd\ + v k a-k, c'_ k = u k a_ k — v k a\ to quasi- 
particle fermion operators a k ,a k , with u k ,v k = [A&± (6 — f cosu;fc)]/(2Afc). The quasiparticle 
energies (12T1) are two-fold degenerate (Afc = A_fc) except for k = n/2 for n odd. 

We can now determine the exact evolution for any n. In the Heisenberg representa- 
tion (dO/dt = i[H,0]), we have a\(t) = e tXkt al(0), a k (t) = e" lXkt a k (0), and the ensuing 
contractions 

(4(*K(t)) = v\ , (a{(t)al k (t)) = -u k v k e 2lXkt , 
with respect to the present initial state (vacuum of the operators c, c') . The average of any 



operator can now be evaluated by substitution and use of Wick's theorem [33 ] . 



1. Evaluation of pit) and C\{t) 

The one-qubit spin flip probability becomes 

p(t) = (c](t) Cj (t)) = Ij^l^k sin2 Xkt , (22) 

n k A k 

where Yl' k = Y^k^i^ 2 denotes integer part). For n = 2, 3 the sum in (1221) reduces 

to a single term [k = 1/2, with uo k = tt/2 and 7r/3 respectively) and we recover exactly Eqs. 
flU and (PJ. 

For n > 4, the evolution of will be in general quasiperiodic. Its upper envelope can 
nevertheless be obtained setting sin 2 X k t = 1 V k in ( 1221) : 

/,n ^ 2 V^' # 2 sin 2 w fc 

p(*) < Pm = - V 77 x 2 , 2 . 2 ' ( 23 ) 

n '—f [b — v cos uj k y + g z sin 
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the maximum of p(t) lying arbitrarily close to p m for sufficiently long time intervals (except 
for rational ratios Xk/Xk')- For low j<d, p m will then exhibit [n/2] peaks, located at 

b = b k = vcosu k , k = §,•■•, [f] -|, (24) 

(i.e. Uk = ir/n, 3ir/n, . . . , (2[n/2] — l)n/n), which are the fields where the quasiparticle 
energies \±k are minimum and vanish for g — > 0. Hence, they are located symmetrically 
around 6 = for even n (o[„/2]-fc = —bk), with a peak at b = (k = n/4) for n/2 odd, but 
asymmetrically for odd n. Moreover, while for b ^ bk, p m oc g 2 , vanishing for g — > 0, at 
b = bk Pm remains finite V (? ^ 0, wrf/i p m ^ 2/n for g — > (Eq. [23]). This implies 

<T - 2^(1 - I) , (25) 

at b = bk for g — > and n > A (and — > 1 for n = 2, 3, 4 as in these cases 2/n > 1/2). Thus, 
by adjusting the field it is always possible to achieve, in principle, finite E\ entanglement 
even for arbitrarily low (but non-zero) values of g. The effect of low anisotropics is just to 
determine the width of these peaks, given by ~ g \ sino; fc | in p m , which increases as g increases 
or as \bk\ decreases. 

The evolution at b = b^ becomes purely harmonic for g — > 0, with 

p(t) | sin 2 X k t , X k = g sin u k . (26) 

The maximum of p(t) is first reached at t k = n/(2g sinuJk), so that the smaller the value of 
g (or u>k), the longer it will take to reach the maximum. In this sense, while the maximum 
entanglement reached in an unbounded time interval is not a continuous function of g for 
g — > at b = bk, that reached in a finite interval [0,tj] will actually vanish for g — * also at 
b = bk, in agreement with the result for g = 0, becoming lower than ( 1251) if t/ < 

The situation at the resonances b = bk is thus similar to that encountered for n = 2 
at b = or for n = 3 at b = v/2. At b = bk the energy gap 2Xk between positive parity 
states with the pair (k, —k) occupied and empty (in particular that between the quasiparticle 
vacuum |0 9 ) and the state o) k a}_ k \Q q )) is minimum, vanishing for g — > (level crossings). Due 
to these degeneracies, at b = bk the aligned state is not approached by any of the Hamiltonian 
eigenstates for g — > 0, remaining distributed over essentially two eigenstates. The previous 
limits (|25|) - (|26|) can then be directly derived from Eq. (120bj) . where for g — > and b = bk, we 
may conserve just the ±k terms in the ^-interaction. The evolution subspace in this limit 
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is then spanned by the original fermionic vacuum |0) (the present initial state) and the two 
particle state \k,—k) = c£c_ fc |0), with ^-independent eigenstates |±) = (|0) =F \k, — k))/y/2 
of perturbed energies ±gsmu>k (i.e., ±Afe). We then obtain (omitting a global phase) 

|^(t)) -> cosA fe t|0) + ism\ k t\k, -k) , (27) 

for the fermionic \^(t)), which leads immediately to Eq. (|26l) . The factor 2/n is just the 
average occupation (cjcj) = J2k'( c k ,c 'k')/ n i n the state \k, —k). 

As g increases, the resolutions of the individual peaks diminish, merging eventually into 
a single broad peak centered at b ph 0. Since the separation between maxima is 5b ph 
(27rv/n) \ sincjfcl, we have the approximate bound g < ttv /n for visible individual peaks. On 
the other hand, it is to be noticed that for n > 5 maximum E\ entanglement can be reached 
only above a certain threshold value g c of g (and then within a certain field window), with 
g c < v V n since at b = and g = v we have exactly p m = (2/n) ^ fc sin 2 u) k = 1/2 for any 
n. In fact, g c ~ v for large n. For (7 ^> (v,b), p m — > 1 (1 — 1/n) for n even (odd), so that 
saturation in C\ is always reached. Finally, for large fields |6| 3> v,g, 

2g 2 s ^' 2 g 2 

k 

implying C™ ~ \f2gj\b\. This asymptotic result is independent of n (for n > 3) and 
coincident with the result previously obtained for n = 3. 

Results for n = 4, 5 and 14, 15 are shown in Figs. [2] and O For n — 4, the resonances 
are located at 6 fc = ±v/y/2, with p m > 1/2 (and hence = 1) for \b 2 — v 2 /2\ < g 2 /2. 
This determines two saturated plateaus in centered at b = b k for g < v , which merge 
into a single plateau centered at b = for g > v . For n = 5 the peaks are located at 
6 fc = v(l ± w 0.81, -0.31, where CJ" -> 2^/5 w 0.98 for p -> (Eq. [25]). Saturation 

is reached only for g/v > 0.67, initially just at the right peak, although for g > v, C™ 
exhibits again a saturated plateau covering 6 = 0. For n = 14 (15), Cj 71 — > 0.7 (0.68) at the 
seven peaks for g — > 0, and saturation is reached for g > 0.92. 

The small or tiny dips in the numerical result for C{" that can be seen in Figs. 2 and 
3 arise due to the occurrence of rational ratios between the quasiparticle energies at 
particular values of b/v, in which case the maximum of p(t) can be lower than the smooth 
upper envelope (|23|) . For instance, for n = 4 the ratio of the two distinct energies X1/2, 
A3/2 becomes 2 at \b\/v = V2{h ± a/16 — 97 2 )/6 (provided 7 < 4/3), where the maximum 
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FIG. 2: (Color online). Maximum concurrence between one qubit and the rest (upper blue curves) 
and between adjacent qubits (lower red+dashed pink curves), reached in the four (left) and five 
(right) qubit chains for two different anisotropies 7 = g/v (type II sectors in C2 depicted again 
with dashed pink lines). For n = 4 the peaks in the global concurrence at b/v = ±l/\/2 are no 
longer resolved for 7 > 1, but remain in the pairwise concurrence. For n = 5, the resonances are 
located at b/v = (l±\/5)/4 and merge again in a saturated maximum for 7 > 1, while the pairwise 
concurrence presents a type II resonance just at the second peak, which again remains visible for 
large 7. Dotted lines in the upper curves depict results obtained with the upper envelope ([23]) . and 
are almost coincident with the numerically obtained maximum in the interval < vt < 40. See 
text for more details. 

reached by p(t) is just (4/5)p m (20% reduction). A reduction in the maximum of p(t) will 
also occur in the vicinity of these values of \b\/v for finite time intervals. This effect gives 
rise to the noticeable dip in C™ at \b\/v « 1.8 for 7 = 1 (the other value \b\/v ~ 0.55 lies 
within the plateau region and its effect on C™ is unobservable) and to those at \b\/v ~ 0.24 
and « 2.12 for 7 = 0.1. 

It should be also mentioned that for short times X^t 1 V/c, p(t) becomes independent of 
n, its series expansion of order m remaining stable for n > m. For instance, up to 0((A&t) 4 ) 
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FIG. 3: (Color online) Maximum concurrence between one qubit and the rest of the chain (upper 
blue lines) and between adjacent qubits (lower red+dashed pink lines) in a n = 14 (left) and n = 15 
(right) qubit chain for different anisotropies, reached in an interval < vt < 180. The dashed pink 
lines depict the maximum of the type II pairwise concurrence, which becomes now lower than the 
type I plateau for 7 > 0.25. Results for C\ obtained with the upper bound ([23]) are also depicted 
(dotted lines, almost overlapping with the blue solid lines). The peaks in C\ are visible for 7 < 0.4, 
and saturation [C\ = 1) is reached for 7 > 0.92. 

in p(t), we obtain, for n > 5, 

pit) w ±g 2 t 2 ll-±t 2 (v 2 + 4b 2 + 3g 2 )], 
d(t) w V2gt[l - ±t 2 {v 2 + Ab 2 + 9s 2 )/24] . 
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It is thus seen that for g ^> (b,g) and n > 8, p(t) exhibits an initial peak at t ~ 1.92/g, 
where p(t) ~ 0.7, with p(t) > 1/2 for 1.2 < gt < 2.75, so that in this limit saturation in C\ 
is rapidly reached (see Fig. The initial peak in C% can be correctly predicted by its 7 th 
order expansion. 

2. Evaluation of C2(t) 

Let us now examine the pairwise concurrence. The relevant elements (1111) of the adjacent 
pair density are 

\ i \, \ /xv 2 sr^'Q 2 cos oa. sin 2 u k . 9 , 
(3{t) = cic J+1 ^ = -p ± *sm 2 A fe t, 

k K 

«(*) = (4(o4 + i(0)o = ^E'^f^ ( 29 ) 

k 

x sin \ k t[ b ~ vc x l SUJk sin A fc t - i cos A fc t] , 

pi(*) = (4(*) c i( t ) c 5+i(*) c i+i(*))o 

= p»(f) -/?*(*) + |o a (f)|, (3()) 

where j < n and in (130]) we have applied Wick's theorem for vacuum expectation values. 

The corresponding results for n = 4, 5 and 14, 15 are also depicted in Figs. [2]l3j It is seen 
that for low g, Czif) presents sharp type II resonances only below the outer peaks of C\, 
and actually just below the rightmost peak for small odd n. In order to understand this 
behavior, we note that for g — > and b = b k , 

(3(t) —> ! cos sin 2 A fc t , \a(t)\ — > ^| sinu k sin2A fc t| . (31) 

These limits can also be directly read from Eq. ( 1271) . as (2/n) cosc^ is the average (<£-Cj-+i) = 
J^ fc/ cosu k '(c k \c k ')/n in the state — fc) whereas a(i) is the average sinu; fc /(c / A |,c / i fc ,)/n 
in the full state ([27]) . The type II maxima of C2 are then obtained for sin 2 \ k t = 1, leading 
to 

C T |[| coscufcl - sin^ fc ^/l - | + ^sin 2 w fc ] , (32) 
in this limit at b = b k . Eq. (1321) is actually positive for 

sin 2 ^<[l-| + ^/(l-^ + ^]-^I + I + 0(^), 
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i.e., u k < 7r/4 or uj k > 37r/4 (\b k \/v > for large n, so that they arise just beneath 

the outer peaks of Ci, the strongest located at the rightmost peak for n odd {k = 1/2) and 
outermost peaks for n even (k = 1/2 or [n/2] — 1/2). Thus, type II resonances in C<i remain 
also finite for g —>■ but are of order nT 1 , becoming smaller than those of C\ for large n 
(C^/Cr oc y/2/n). The scaled concurrence nC™ remains nevertheless finite for large n. 

For n = 3 we exactly recover from fl3"2"l) the previous result = 2/3 for the type II 
peak. For n = 4, Eq. fl32D yields = (2y/2 -l)/4a 0.46, whereas for n = 5 it leads to 
a single peak at u k = 7r/5, of height m 0.41. For n = 14, there are sharp type II peaks at 
the outer resonances, of height « 0.22, plus smaller peaks at the next resonance, of height 
~ 0.08, which rapidly fall below the type I plateau. For n = 15 the visible type II peaks are 
asymmetric and appear at b k /v m 0.98,0.81 and —0.91, with heights ~ 0.21,0.08 and 0.15. 

For g — > there are also type I maxima of Ci at b = b k , visible in the central region (Fig. 

|3]). These maxima are broader and occur at times determined by 

1 — - sin 2 u)u 
cos(2A fc t) = , 

ysin 2 uj k + (1 - | sin 2 u k ) 2 
(the first peak at t\ ~ ir/(%\ k ) for u k ~ 7r/2), where the concurrence approaches for g —>■ 
the value 



- iiy/0- - ~ sm 2 ^,)2 + sin 2 ^ - (1 - ^sin 2 u; fc )] . (33) 

Since this is an increasing function of | smu>k\, i.e., a decreasing function of \b k \, the type / 
maxima fall below those of type II for low | sinc^l (| sincijfc| < 0.66 or \b k \/v > 0.75 for large 
n). Moreover, at the highest type I peak (tu k ~ vr/2), C™ ~ 2(v2 — l)/n for large n, which 
is just 21% of the highest type II peak (C™ ~ 4/n). For n = 3 we also recover from ( 133]) the 
previous exact result C™ = 1/3 in the type I plateau, while for n — 2 it yields the correct 
maximum value C™ = 1. For n = 4 and 5 we obtain ~ 0.14 and C™ ~ 0.07, 0.2 at the 
type I peaks, while for n = 14, 15, ~ 0.07, 0.06 at the centermost type I peak for g — > 0. 

As g increases, the lower type I resonances in C™ become rapidly smoothed out, merging 
into a broad plateau (Figs. |2|3|) . Moreover, while for low n the type II peaks remain visible 
even for large g (Fig. EJ), as n increases these peaks become as well superseded by the type 
I plateau (Fig. [3]), which is discussed below. On the other hand, for \b\ 3> v,g, we obtain, 
up to first order in g/\b\, v/\b\, C 2 (t) w 2|a(t)| < C7 2 m , with 

n \ b \ k \ b \ 

16 








1.0 
U 
0.5 I- 



-i — i — i — | — i — i — i — | — i — i — i — r 
7 = 0.1 b/v = 2 



i 



-1 I I I I I I I I I I r 
7=1 b/v = 2 




10 20 30 
gt 



10 20 30 
gt 



FIG. 4: (Color online). The evolution of C±(t) (upper curves in blue) and C2(t) (lower curves, in 
red and pink) for n = 15 at two different anisotropies and different fields. The central panels depict 
the evolution at the outer resonance b c /v = cos(-7r/n) « 0.98. Both the type I (red) and type II 
(pink, dashed lines) sectors of C^if) are indicated. 

in agreement with the previous result for n = 3. In this limit, « 



3. Temporal Evolution 

Fig. H] depicts C\(t) and C2(t) for n — 15 at two different anisotropies, at and away from 
resonances. For low 7 (left panels), we observe a low frequency periodic-like evolution of 
Ci(t) and Czit) at the outer resonance (b/v ~ 0.98), in agreement with ([26]) and (I3"T1) . with 
Czit) exhibiting regions of both type I and type II entanglement, whereas for large fields 
b = 2v both Ci(t) and C^t) become very small, with C^t) of type I. Both C\{t) and C*2(t) 
are also smaller for 6 = (with C2 (t) again of type I) , which here corresponds approximately 
to a minimum of CJ" and C™. 
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5 10 
gt 

FIG. 5: (Color online). Evolution for large anisotropy and short times of C\(t) (upper curves in 
blue) and C2(t) (lower curves, in red), for neighboring odd-even systems. 

On the other hand, for 7 = 1 the emerging global entanglement is non-negligible for 
all moderate fields, with saturation in C\ reached for b < v. In this case C^if) does not 
follow the behavior of C\{t) for low fields, where it strictly vanishes at finite time intervals, 
although for large b > v the evolution of C^t) becomes again similar to that of C\{t) (with 
C™ ~ C™/y/2), and intervals of vanishing value are removed. Thus, the average pairwise 
entanglement is in this case enhanced by a large field b pa 2v, in comparison with that for 
b pa v , as a consequence of the lower global entanglement. In other words, the decoherence 
of the pair for large 7 due to the interaction with the spin chain (representing here the 
environment for the pair) is prevented by large fields. 

It is also seen that the evolution for 7 = 1 (g = v ) and b = is strictly periodic. In this 
case Afc = v V k and Eqs. ( 1291) become independent of n for n > 4 and of the form 

p(t) = I sin 2 v t , (3{t) — 0, a(t) — — ij sin2wt , 
Ci(t) = I sinvt\y/2- sin 2 vt , (34) 
C 2 (t) = |sint)t|Max[|cost)t| - | sinvt|/2, 0] . (35) 

Hence, C\{t) reaches saturation when | sin ft | = 1, whereas has maxima when cos 2vt = 
l/y/E, where (^(t) = (y/E — l)/4« 0.31, and vanishes in the interval where | cos Vt\ < 1/y/E 
or when sin ft = 0. The previous maximum of C2 is already close to the maximum obtained 
for large 7 (see below) and is higher than the resonant values for n > 9. 

Fig. dSJ) depicts the typical evolution for short times and large anisotropy. As seen here, 
the plateau in the maximum concurrence C™ arising for g > (v, b) is originated by the first 
maximum in the evolution of Ci(t), which exhibits in this region a prominent initial "burst" 

18 



followed by intervals of vanishing value (i.e., decoherence of the pair) and lower revivals 
(near the most prominent minima of C\(t)). For g ^> (b, v) and n > 5, the initial peak 
of C*2 occurs at gt ~ 0.66, with height C™ ~ 0.35, and is practically independent of n. 
The resonances in C™, of order n _1 , become then rapidly covered by the plateau as n or g 
increases. This initial peak can be approximately reproduced by a fourth order expansion 
of C2(t), given for n > 5 by 

C 2 (t) « - fpt - I* V + & 2 + 3<? 2 ) 
+^t 3 (26 2 + 3/-t; 2 ). 

Nonetheless, odd-even differences and n-dependence do arise for longer times (gt > 10 in 
the case of Fig. [5]) and affect the revivals of C^- 

Let us finally mention that as the resonances arising for low 7 develop their first maximum 
at tk = 7r/(2<7 sin 00%), the relevant timescale for their observation is r ~ h/(^v) « 17/7, where 
t v ~ H/v is the operation time associated with the hopping strength w, and should be smaller 
than the characteristic decoherence time of the chain determined by its interaction with 
the environment. This limits the smallness of the anisotropy (i.e., 7 > t v /t<i) and hence 
the sharpness of the peaks. For instance, if 7 = 0.1 and v ~ 0.02 meV, which is a typical 
strength for realizations based on quantum dots electron spins coupled through a cavity 



mode 171 ]. r~3x 10 _10 s, which is smaller than the typical decoherence time 17|. On the 
other hand, the results for C2 represent the evolution of the entanglement of an adjacent 
pair in the present spin chain environment, and indicate that resonances remain finite at 
the pairwise level in such scenario. 



IV. CONCLUSIONS 



We have examined the entangling capabilities of a finite anisotropic XY chain with 
constant parameters for an initially completely aligned state in the transverse direction. 
The exact analytical results obtained (valid for all n) show that the maximum attainable 
entanglement exhibits for low anisotropy 7 a clear resonant behavior as a function of the 
transverse magnetic field, with peaks at those fields where the effective quasiparticle energies 
Afc are minimum and vanish for 7 = 0. At these fields, the energy levels become then 
degenerate for 7^0 and the aligned state remains mixed with its degenerate partner for 
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arbitrarily small but non-zero 7. The height of these resonances remains thus finite for 
7^0 and their width is proportional to the anisotropy, implying a fine field sensitivity apt 
for efficient control, although the time required to reach the peak is proportional to 7 -1 and 
the height decreases as the number of qubits increases. The resonances are notorious in the 
maximum global entanglement between one-qubit and the rest of the chain, and are present 
as well in the entanglement of other global partitions. 

They also arise in the maximum pairwise concurrence, and can be of both spin parities, 
although they are of lower height and decrease more rapidly with n, being hence more 
easily smoothed out for increasing 7. Here we have shown that type II (I) resonances 
become dominant at large (low) critical fields for adjacent pairs, those of type II being 
extremely narrow. Another feature is that odd-even differences in the resonant behavior 
remain appreciable for moderate n, odd chains exhibiting field sign sensitivity both in the 
global and pairwise peaks. On the other hand, saturation can be reached in the global Ei 
entanglement within a certain field window above a threshold anisotropy (7 ~ 1 for large 
n), but not in the pairwise entanglement, whose maximum exhibits instead a broad low 
plateau for large 7 and hence low field sensitivity. Let us finally remark that resonances 
of the present type will also occur for non-adjacent pairs as well as for other geometries or 
interaction ranges, although details (i.e., relative widths and strengths) may certainly differ 
from the present ones and are currently under investigation. 

RR acknowledges support of CIC of Argentina. 
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